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'~~^ • Abstract 

OO 

CN| I We present an unexpected connection between two map enumeration prob- 

lems. The first one consists in counting planar maps with a boundary of 

Q ■ prescribed length. The second one consists in counting planar maps with two 

^ . points at a prescribed distance. We show that, in the general class of maps with 

r^ ! controlled face degrees, the solution for both problems is actually encoded into 

'^ [ the same quantity, respectively via its power series expansion and its continued 

fraction expansion. We then use known techniques for tackling the first prob- 
lem in order to solve the second. This novel viewpoint provides a constructive 

(^ . approach for computing the so-called distance-dependent two-point function 

^ i of general planar maps. We prove and extend some previously predicted exact 

formulas, which we identify in terms of particular Schur functions. 
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1. Introduction 

1.1. General introduction 



Maps are fundamental objects of combinatorial theory, introduced by Tutte in the 
^ ' 60's 0,11, which also appeared later as natural models for random surfaces in physics 0. 

j^ ■ Many questions about maps boil down to enumeration problems, which in turn received 

a lot of attention from different communities. Different techniques of enumeration have 
been developed, for instance combinatorialists' quadratic method [Q, physicists' matrix 
integrals and more recently bijective coding by trees |]§-|^. Most results deal with 
global enumeration problems, say counting families of maps with a control on their size 
and topology. The coding by trees allowed to address more refined questions about 
the distance in maps, say counting maps with a prescribed radius []T^ or with marked 



pointed at prescribed distances |Tl[]. This approach led to important applications in 
probability theory for the rigorous construction of scaling limits of large random maps 

One of the prominent global quantities, considered both in the recursive decom- 
position and the matrix integral approach, is the number of maps with one boundary, 
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with a control on both the total map size and the boundary length. More precisely, 
this information is captured by the moments, which are generating functions of maps 
with a fixed boundary length. These moments are determined by a closed system of 
equations, called either Tutte's or loop equations respectively in the combinatorial or 
physical communities. 

On the other hand, one of the simplest observables related to distance is the so- 
called two-point function, which is the number of maps with given size and having 
two marked vertices at fixed distance. In particular the dependence of this two-point 
function over the distance gives directly the average profile of random maps. The 
distance-dependent two-point function was first considered in [|T3[ in the case of planar 



triangulations, where its scaling form was predicted via a "transfer- matrix" approach in 
the spirit of Tutte's decomposition. An exact discrete expression for the two-point func- 
tion was given in [|TT| for planar quadr angulations and more generally bipartite planar 
maps with a control on face degrees, whose scaling form agrees with |jT^. This exact 



expression makes use of a coding of such maps by trees, leading to discrete recurrence 
equations whose solution was guessed. 

The purpose of this paper is to present an unexpected yet remarkable connection 
between the two above notions showing that the information about the distance in maps 
is actually hidden in the global problem of enumerating maps with a boundary. More 
precisely, the moments and the distance-dependent two-point function constitute two 
possible expansions of the same quantity, the resolvent: the moments form its power 
series expansion while the two-point function encodes its continued fraction expansion. 
Using the standard theory of continued fractions this allows in particular to obtain 
explicit expressions for the two-point function via the known techniques for computing 
the moments. This program is carried out in detail in this paper. We begin by an 
overview of our main results in Sect. 1.2, where precise definitions for the moments and 
two-point function are given. The actual connection between them is established in 
Sect. 2 by use of the coding of maps by appropriate trees called mobiles. In Sect. 3 
we derive an explicit expression for the moments in terms of weighted paths by various 
techniques: a direct non-constructive check, the solution of Tutte's equation and a 
bijective construction. This third approach also gives rise to a one-parameter family 
of expressions for the moments forming so-called conserved quantities. Our expression 
for moments is used in Sect. 4 to derive explicit expressions for the two-point function 
in terms of Schur functions. The special case of bipartite maps, which was the scope 
of [|11|, is discussed in Sect. 5, while the particularly simple cases of triangulations 



and quadr angulations are addressed in Sect. 6. Concluding remarks and discussion are 
gathered in Sect. 7. 

1.2. Overview of the main results 

A planar map is a connected graph (possibly with loops and multiple edges) drawn 
on the sphere without edge crossings, and considered up to continuous deformation. It 
is made of vertices, edges and faces. The degree of a vertex or face is the number of 
edges incident to it (counted with multiplicity) . A map is rooted if one of its edges is 



distinguished and oriented, the root face being the face on the right of the root edge 
and the root degree being the degree of the root face. 

In this paper, we generaUy consider enumeration problems for rooted planar maps 
subject to a control on face degrees, i.e. for each positive integer k we prescribe the 
number of faces with degree k. In the language of generating functions, this is equiva- 
lent to attaching a weight gk to each face of degree /c, where {gk)k>i is a sequence of 
indeterminates, so that the global weight of a map is the product of the weights of its 
faces. Setting gk = for all odd k amounts to considering the planar maps which only 
have faces of even degree, which are the bipartite planar maps. Drastic simplifications 
occur in the bipartite case and these will be discussed in due course. 

As a general preliminary remark, most quantities introduced in this paper will 
be formal power series in the gk^s (with integer or rational coefficients) but, for the 
sake of concision, this will not be apparent on the notations - we shall write X rather 
than X{gi, (72, (73, . . .) for instance. If extra variables are involved, they shall be written 
explicitly. 

The fundamental observation of this paper is a combinatorial identity between 
some a priori unrelated families of generating functions for rooted planar maps with the 
above weights. On the one hand, for each positive integer n, we consider the generating 
function Fn for rooted planar maps with root degree n, which we call the nth moment. 
By convention we do not attach a weight g^ to the root face and we set Fq = 1. We 
may combine all values of n into the generating function F{z) = Xl^o ^n^^ where 
z is an extra variable. This quantity essentially coincides with the planar "resolvent" 
encountered in the context of matrix integrals. On the other hand, again for each 
positive integer n, we consider the generating function for planar maps with two marked 
vertices at distance n, which we call the distance-dependent two-point function. To 
properly root these planar maps, we take as root an edge incident to and pointing 
away from one of the marked vertices. Then the endpoint of the root edge must be at 
distance n — 1, n or n + 1 from the other marked vertex. We distinguish between these 
three cases and denote by r^, t^ and r^+i the respective generating function (clearly 
reversing the orientation of the root edge allows to identify the third case with the first 
one for n -^ n + 1). We furthermore introduce the "cumulative" generating functions 
Rn = X]r=i ^* ^^^ -^^ ~ Sr=o^«' corresponding to two marked vertices at distance 
less than or equal to n. For convenience, a conventional term 1 is also included in ri. 
Then, our observation is that the sequences {Rn)n>i and {Sn)n>Q = {VT^)n>o form 
the Jacobi type continued fraction (or J-fraction) expansion of F{z), namely 

00 ^ 

F{z) = J^FnZ- = — , . (1.1) 

1-^0. ""'' ^ , 

1-S,z- ^^^ 
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As discussed in detail in Sect. 2, this follows from the correspondence between pla- 
nar maps and some labeled trees called mobiles H, and the combinatorial theory of 



continued fractions |1I4| , ^5(] . Calling Fn the nth moment is consistent with the usual 



terminology of J-fractions |T^. A classical result states that the sequence {Rn)n>i is 



closely related to the Hankel determinants of moments, namely 



-n„_i 0<i,j<n 



n- 



(with the convention H-i = 1) and similarly the sequence {Sn)n>o is given by Hankel 
minors through 

^ H ~ 

y^S^ = -Fr^ ^ri = det Fi+j+5^ (1.3) 

^—Z rtn 0<i,j<n 

1=0 

In the specific context of planar maps, we derive in Sect. 3 a general formula for the 
moments F^ in terms of the weights {gk)k>i- A very peculiar structure emerges when 
this formula is substituted into the determinants H^ and Hn' after some simple ma- 
nipulations we recognize instances of a classical identity for symplectic Schur functions, 
see Sect. 4 for details. Let us now state the important results of Sects. 3 and 4. 

Our expression for the moments involves generating functions for three-step paths: 
these are lattice paths in the discrete Cartesian plane which consist of up-steps (1, 1), 
level-steps (1, 0) and down-steps (1, —1). We denote by P{n; R, S) the generating func- 
tion for three-step paths going from (0, 0) to (n, 0), where a weight 5' is attached to each 
level-step, a weight \/R is attached to each up- or down-step and the global weight of a 
path is the product of the weights of its steps (clearly, this global weight involves an inte- 
ger power of R). Furthermore, we say that a path is positive if it only visits vertices with 
non- negative ordinate, and we denote by P~^{n; i?, S) the generating function for positive 
three-steps paths from (0, 0) to (n, 0), also known as Motzkin paths of length n, with the 

same weights. We have the explicit expressions P(n; R, S) = ^|"o (i\)'^(n-2i)\ ^^ ^^~'^^ 

and P+{n;R,S) = E1=o^^ W+iW^^2ly.^'^''~^'- Then, the general formula for the 
moments is 

oo 

F^ = Y,A,P+{n + q;R,S) (1.4) 

where the coefficients Ag are given by 



A, = Ri 6,,o - Yl 9kP{k -q-2;R,S)\ (1.5) 

\ k=q+2 J 

and where the step weights R, S are themselves power series in {gk)k>i, which are 
implicitly determined by 

oo 1 °o c2 

S = J2 9kP{k - 1; i?, ^), ^ = 1 + 2 E 3kP{k; R,S)-—. (1.6) 

fe=l k=l 



Now, substituting the expression ( |1.4| ) into the Hankel determinant H^, elementary 
row and column manipulations, which amount to a natural decomposition of lattice 
paths, lead to 

oo 

n{n-\-l) ^ — ^ 

Hr, = R—^ det (S,_, -S,+,+2), B, = y^A^P,{q:R.S) (1.7) 

where Pi{q;R,S) is the generating function for three-step paths from (0,0) to {q,i). 
Note that P-i{q; -R, S) = Pi{q; -R, 5"), hence B-i — Bi. Similarly for H^ we have 

~ Ti'^ + n + 1 

Hr, = in + l)SHn + R—^ det (S.-,-5,.„ - S,+,+5^,„+2). (1.8) 

At this stage, we fix a positive integer p and consider maps whose faces have degree at 
most p + 2, i.e. we set gr/j = for /c > p + 2. This implies, by the expression ( [1.5|) , that 
Aq vanishes for q > p and moreover, since Pi{q; R, S) = for |z| > q, that Bi vanishes 
for i > p. The matrices appearing in (|1.7|) and ( p..8|) are therefore band matrices. Their 
determinants can be identified, up to a normalization, with characters of the symplectic 
group Sp2p [0 , also known as symplectic Schur functions. Following |T^ , we denote 



by sp2p(A,x) the symplectic Schur function associated with the partition A, which is 
a symmetric function of 2p variables x = (xi, 1/xi, X2, 1/^2, . . . , Xp, l/xp). Here, these 
variables are the 2p roots of the "characteristic equation" 



p 



Br.x'' = (1.9) 



-p 



(recall that -B-n = B^, hence the equation is invariant under x — )■ 1/x). By the 
expressions (|1.5| ) for Aq, ( pT7| ) for Bi and the identity X]i=-q PiiQ'i -^5 S) ^* = {VRx + 
S + y/R/xY, the characteristic equation reads explicitly 

p+2 k-2 / /^^ 



l = 5^(7fc5^P(A;-2-Q;i?,5) v^x + ^ + ^— . (1.10) 



X 
k=2 g=0 



The partitions associated with the determinants appearing in (|1.7| ) and ( |1.8| ) are re 



spectively the "rectangular" partition Ap,n+i — {n + 1)^ made of p parts of size n + 1, 
and the "nearly-rectangular" partition Ap^^+i = (n + lY~^n made of p — 1 parts of 
size n -(- 1 plus one part of size n. In summary, we find that Hn oc sp2p(Ap,n+i, x) and 
Hn — {n + l)SHn oc sp2p(Ap^„_|_i, x). Using the Weyl character formula for sp2p(A, x), we 
end up with the following "nice" formulas for R^ and Sn, which involve determinants 
of size p, independently of n: 

det ix1+^+' - x-""^-') det {xr^+' -x-''+^-') 
R^ = r'-^^^^ '-^^^^ -^ , (1.11) 

det (x^ ■' — x~^~' 



det (x:+'+^'-'-^-a;p-'-^'+'-^) det (a;:+^'-'-^ - x-"-^'+'-^) \ 

s^ = s-Vr I '-^^^^ -^. 3^ '-^^^^ -. ^ 

, det (x" ^-x. " ^) det (x^ ^ -X. " •^) 

(1.12) 
A number of simplifications occur in the bipartite case, discussed in detail in Sect. 5. 
Clearly F^ = for n odd and the continued fraction expansion of F{z) naturally reduces 
to the Stieltjes type 

F{z) = Y, F2nz''' = ^^. (1.13) 



n=0 1 



Riz' 



1 

Consistently, the quantities 5", H^ and Sn vanish while the Hankel determinants factorize 
as 

Hr, = hfl,h[l,,, /ii°)= det F2.+2,, /ii^)= det F2,+2,.+2. (1.14) 

L2"J L~2~J 0<i,j<n 0<i,j<n 

This factorization property is also apparent at the level of Schur functions. Indeed, for 
bipartite maps whose faces have degree at most 2p + 2, we have a characteristic equation 
of the form 



Eb= 



i2,.i""=0 (1.15) 

n^—p 

whose roots are of the form (xi, 1/xi, — xi, — 1/xi, . . . ,Xp, 1/xp, — Xp, — 1/xp), and 
sp4p(A, x) can be factorized as the product of two symmetric functions of the vari- 
ables (xf, 1/x^, . . .Xp, 1/Xp) (a symplectic and an odd orthogonal Schur function). In 
the end, Rn can be expressed in terms of determinants of size p rather than 2p, which 



matches the expression found in ill 



The cases of triangulations and quadrangulations are particularly simple as 
the roots of their characteristic equation are of the respective forms (x, 1/x) and 
(x, 1/x, — X, — 1/x). Consequently, the two-point functions in these cases can be ex- 
pressed in terms of a single quantity x solution of some algebraic equation. This is 
discussed in Sect. 6 and a simple interpretation via one-dimensional hard dimers is 
given. 

2. Prom planar maps to continued fractions, via mobiles 

The main purpose of this section is to establish that the generating functions for 
planar maps F^, Rn and Sn defined in Sect. 1 are related by the J- fraction expansion 
( p. . 1|) . The plan is as follows. We recall some general results of the combinatorial theory 
of continued fractions (Sect. 2.1), before reviewing the coding of planar maps by mobiles 
(Sect. 2.2). We then interpret Fn, Rn and Sn as generating functions for mobiles and 
complete our proof (Sect. 2.3). Finally we mention a few other outcomes of our approach 
(Sect. 2.4). 
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2.1. Reminders of the combinatorial theory of continued fractions 

In his seminal paper fl^ , Flajolet gave a combinatorial interpretation of continued 



fractions in terms of Motzkin paths. Let us recall a few important results of this the- 
ory. For the purposes of this subsection, {Rm)m>i and {Sm)m>o may denote arbitrary 
sequences of elements of a commutative ring. That is to say, we may temporarily forget 
about their map-related definition given in Sect. 1. We consider Motzkin paths, and 
more generally positive three-step paths as defined above, with the following ordinate- 
dependent weights: for all m, 

- each down-step of the form (t, m) — t- (t + 1, tti — 1) receives a weight Rm {'m > 1) 

- each level-step of the form (t, m) — t- (t + 1, m) receives a weight Sm (?^ > 0) 

- each up-step receives a weight 1. 

Note that attaching ordinate-dependent weights to up-steps would add essentially no 
generality. Let us now define F^ as the generating function for Motzkin paths of length 
n with ordinate-dependent weights (we also temporarily forget about the map-related 
definition of Fn) and, attaching a further weight z per step (of any type), F{z) = 
X]^o ^nZ"' the generating function for Motzkin paths of arbitrary length. Then, the 
Continued Fraction Theorem (as named in |T^) states that F{z) is given by the J- 
fraction (|1.1|). 



In order to recall of its derivation, let us introduce a few notations which also prove 
useful later. Given arbitrary non-negative integers d, d' and n, we denote by Zd^d'i'n) 
the generating function for positive three-step paths from (0, d) to {n^d'), with the 
above ordinate-dependent weights. Furthermore, we denote by Z'^^,{n) the generating 
function for such paths restricted to have all their ordinates larger than or equal to 
min((i, d') (instead of 0). Clearly F^ = ZQ^Q{n) = ZQQ{n). Now, for the case d = d' we 
have the relation 

^ ' l-SdZ- Rd+1 z^ E ZJ^i,d+i (n) ^" 

~ n>0 

which simply translates the "arch decomposition" , namely by cutting a path contribut- 
ing to the left-hand side at each occurrence of the ordinate d, it is bijectively decomposed 
into a sequence of two types of objects: either (i) level-steps at ordinate d, weighted by 
Sd, or (ii) arches made of the concatenation of an up-step, a restricted path starting and 
ending at ordinate d+1 and a down-step, with an overall weight Rd+iZd+i^d+iin') for 
some n' . By iterating (|2.1|) starting at (i = 0, we deduce (|l.lj). Note that the J- fraction 



is a well-defined power series in z, because computing F^ only requires Ln/2J iterations 
(as Motzkin paths of length n reach at most the ordinate [n/2\ ) . If we iterate ( |2.1|) 
starting at an arbitrary d, we obtain 

E^J»^" = ^-^-—2 (2-2) 



'->o 1 - Sdz 



Rd+2Z 

Sd+iz - — — 



where the right-hand side is caUed a truncation of the fundamental fraction ( |1 . 1|) . Many 
other identities are known (e.g for X]n>o ^dd' ('^) ^^ ^^^ X]n>o ^d,d' (n) z^ via "last- 
passages decompositions" ) but we shall not need them here. 

A particular case is when the weights do not depend on the ordinate, i.e. we set 
Rm — R for all m > 1 and Sm = S for all m > 0. Then, we recover the three-step path 
generating functions defined in Sect. 1, namely Zd,d'in) = R^'^~'^ )/^Pd'-d(^; -R, S) and 



'd,d' 



(n 



Rid-d')/2 p+ ^^^^. jl^ gy Note that the i?(<^-^')/2 prefactor arises from the 



slightly different weighting convention. Most notably, we have 



J2P^in;R.S) 



n>0 



1-Sz 



Rz' 



(2.3) 



1-Sz 



Rz' 



1 



Clearly this series is the solution of the quadratic equation X = 1 + SzX -\- Rz^X^. 
More generally, if we only assume that Rm — ^ R and 5'^^ — )■ 5" for rri — ;■ oo (in any sense 
of convergence), then obviously Zd^di'n) — )■ P{n;R,S) and Z^^{n) — )■ P~^{n;R,S) for 
(i — i- oo. 

Let us now consider the inverse problem of determining the sequences (-Rm)m>i 
and {Sm)m>o in the J-fraction expansion, knowing the power series expansion of F{z) 
i.e. the sequence {Fn)n>o- It is a classical result that this problem is solved using Hankel 
determinants, and the solution is given by relations ( |1.2| ) and ( [1.3| ). See for instance 
the beautiful proof by Viennot |T^ based on a combinatorial interpretation involving 
configurations of non-intersecting Motzkin paths. 



2.2. Review of the coding of planar maps by mobiles 



y\ 



(a) 



@ 

• 4 • 

® 

(b) 



Fig. 1: The mobile construction applied (a) to an edge of type (n, n — 1) 
and (b) to an edge of type (n, n). In (a) a mobile edge (thick line) connects 
an unlabeled vertex (big dot) to a labeled one. In (b) a bivalent flagged 
vertex (lozenge) is created and connected to both adjacent unlabeled ver- 
tices. 
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We now return to the context of planar maps and review their coding by mobiles 
presented in Ref. ^. The coding applies to pointed uisips, i.e. maps with a distinguished 
vertex called the origin, and works as follows: we start by labeling each vertex of the 
map by a non-negative integer equal to its graph distance from the origin and add a 
new unlabeled vertex at the center of each face. With the above labeling, any edge of 
the map is either of type (n, ti — 1), i.e. connects a vertex labeled n to a vertex labeled 
71 — 1 for some n > 1, or of type (n, n) for some n > 0. For each edge of type (n, ti — 1), 
we draw a new edge connecting its incident vertex labeled n to the unlabeled vertex 
sitting at the center of the face on the right of the edge (oriented from n to n — 1). This 
procedure is displayed in Fig. l-(a). For each edge of type (n, n), we add in its middle 
a "flagged" vertex with flag n and connect it by two new edges to the two unlabeled 
vertices at the center of the two incident faces, as displayed in Fig. l-(b). It was shown 
in IP] that the new edges form a tree which spans all the labeled vertices with label 
n > 1 (i.e. all the original vertices of the map but its origin), as well as all the added 
unlabeled vertices and all the added flagged vertices. 



5 
4 
3 





Fig. 2: An illustration of the property (P) around an unlabeled vertex. 



The labels and flags satisfy the property: 

(P) for each unlabeled vertex, the clockwise sequence of labels and flags on its adjacent 

vertices matches exactly the sequence of labels and flags obtained by attaching a 

label n to each down-step (t, n) — )■ (t -|- l,n — 1) and a flag n to each level-step 

(t, n) — )• (t -|- 1, n) of some three-step path with identical initial and flnal ordinates. 

Note that the three-step path in (P) is unique up to cyclic shifts, and that its length is 

given by the number of flags plus twice the number of labels. Returning to the original 

map, the sequence of ordinates in the three-step path simply reproduces the clockwise 

sequence of distances from the origin to the vertices around the face associated with the 

unlabeled vertex at hand, while the length of the path is nothing but the degree of this 

face (see Fig. 2). 



In all generality, a mobile is defined as a plane tree with three types of vertices: 
unlabeled vertices, labeled vertices carrying an arbitrary integer label and flagged ver- 
tices carrying an arbitrary integer flag, with edges connecting only unlabeled vertices 
to labeled or flagged ones in such a way that the flagged vertices have degree 2 and 
property (P) holds around each unlabeled vertex. Note that shifting all labels and flags 
in a mobile by a flxed integer preserves property (P) so that the resulting object is still 
a mobile. 

The above coding associates to each pointed map a mobile. As shown in P], 
it is a bijection between pointed planar maps and mobiles satisfying the two extra 
requirements: 

(Rl) labels and flags are respectively positive and non-negative; 
(R2) there is at least one label 1 or one flag 0. 
Strictly speaking, the tree reduced to a single vertex labeled 1 satisfles (Rl) and (R2), 
and corresponds to a conventional empty map. 

Furthermore, the coding yields the following one-to-one correspondences between 
elements of a map and of its associated mobile. Map vertices at distance n > 1 from the 
origin in the map correspond to mobile vertices labeled n. Map edges of type (n, n), i.e. 
connecting two vertices at distance n from the origin, correspond to flagged vertices with 
flag n (n > 0) while map edges of type (n, n — 1) correspond to mobile edges incident to 
a labeled vertex with label n (n > 1). Faces of the map correspond to unlabeled vertices 
of the mobile and moreover, the clockwise sequence of distances from the origin of the 
vertices incident to a given face is directly read off the ordinates of the (unique cyclic) 
three-step path ensuring property (P) around the corresponding unlabeled vertex. 

The mobile coding a pointed map is unrooted, i.e. it has no distinguished edge 
nor vertex. We deflne conventionally a rooted mobile as a mobile with a distinguished 
edge incident to a labeled vertex (the other incident vertex being necessarily unlabeled) , 
whose label is then called the root label. The mobile reduced to a single labeled vertex 
is considered as a rooted mobile even if it has no edge. Regarding flagged vertices, we 
flnd more convenient to introduce the notion of half-mobile, which is deflned exactly as 
a mobile except that it has one particular flagged vertex of degree 1 (the other flagged 
vertices being of degree 2 as before) whose flag is called the root flag. A mobile having a 
distinguished edge incident to a flagged vertex with flag n can clearly be seen as a pair of 
half-mobiles having the same root flag n. Rooted mobiles and pairs of half-mobiles code 
naturally pointed rooted maps, namely maps with both a marked vertex (the origin) 
and a marked oriented edge (the root edge). More precisely, a pointed rooted map of 
type (n — )■ n — 1), i.e. whose root edge points from a vertex at distance n from the 
origin to a vertex at distance n — 1, is coded bijectively by a rooted mobile satisfying 
(R1)-(R2) and having root label n (we distinguish the mobile edge associated with the 
root edge), while a pointed rooted map of type (n — t- n) is coded bijectively by a pair 
of half-mobiles with root flag n and whose union also satisfles (R1)-(R2) (we split the 
mobile at the flagged vertex associated with the root edge, whose orientation allows to 
distinguish the two half-mobiles). For the last remaining type (n — 1 — )■ n) we may 
simply reverse the root edge orientation. This encompasses all possible types of the 
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root edge, when n varies. Note that rooted maps can be seen as pointed rooted maps 
of type (1 ^0) or (0^0). 

It is worth noting that a bijection also exists with mobiles where the conditions 
(Rl) and (R2) are waived. It is obtained by now considering pointed rooted maps where 
the distance between the origin of the map and the starting point of the root edge is 
no longer fixed. Denoting by d the distance of the extremity of the root edge farthest 
from the origin, we now label each vertex by its distance from the origin minus d so 
that the pointed rooted map is either of type (0 — )■ 0), (0 — )■ —1) or (—1 — )■ 0). Applying 
the mobile construction rules of Fig. 1 to all edges of the map creates a mobile where 
the condition (P) still holds but where the conditions (Rl) and (R2) have been waived. 
For type (0 — )■ 0) we finally obtain a pair of half-mobiles with root flag 0, while for 
type (0 —7- —1) or (—1 — )■ 0) we obtain a rooted mobile with root label (for a proper 
bijection we need to adjoin a sign ±1 to the rooted mobile, in order to keep track of the 
orientation of the root edge). 

2.3. Mobile interpretation of the generating functions F^, Rn, Sn 

We now want to show that ( |lT|) holds, when F„, Rn and Sn are the generating 
functions for planar maps defined in Sect. 1. This requires identifying these quantities 
as generating functions for mobiles, and Fn will naturally appear as a sum over Motzkin 
paths with weights Rm per down-step starting at ordinate m and Sm per level-step at 
ordinate m. As a general preliminary remark, note that the weight g^ per face of degree 
k in the map translates, in the mobile language, into a weight gk per unlabeled vertex 
whose three-step path associated via (P) has length k. 

To begin with, let us discuss the case of Rn {n > 1). From our definition of 
Sect. 1.2 and the above discussion, Rn = J27=i ^« enumerates pointed rooted maps of 
type (z — 7- z — 1) for some i between 1 and n, which are in bijection with rooted mobiles 
satisfying (R1)-(R2) and having a root label i between 1 and n. These mobiles are 
in turn in bijection with rooted mobiles satisfying (Rl) only and having root label n. 
Indeed, we may transform the label i in the first family of mobiles into a label n by 
shifting all labels and fiags of the mobile by n — i. The condition (Rl) is preserved by 
this non-negative shift but (R2) is no longer valid if z < n. Conversely, assuming that 
(Rl) is satisfied in a rooted mobile with root label n > 1, the minimal label among 
labeled vertices necessarily lies between 1 and n, i.e. is of the form n — i + 1 for some 
i between 1 and n. From the general condition (P) around each unlabeled vertex, we 
easily deduce that the minimal fiag in the mobile is at least n — i. Shifting all labels 
and fiags by the (non-positive) quantity i — n preserves (Rl) and restores (R2), while 
it transforms the root label n into an i between 1 and n. To summarize, Rn may be 
identified as the generating function for rooted mobiles satisfying (Rl) only and having 
root label n. In this identification, the conventional term 1 added in ri, hence in Rn, 
accounts for the mobile reduced to a single labeled vertex with label n. 

As for Tn = ^iLo^J' i^ enumerates pointed rooted maps of type (z — )■ i) for some 
i between and n, which are in bijection with pairs of half-mobiles with root fiag i 
and whose union satisfies (R1)-(R2). By a similar argument as above, these are in 
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bijection with pairs of half-mobiles with root flag n satisfying (Rl) only. Note that (Rl) 
independently applies on each half-mobile. We arrive at T^ = (3^)'^ where we identify 
directly S^ as the generating function for half-mobiles satisfying (Rl) only and with 
root flag n. 



m-v 




external face of degree n 




(«,0) 



ZcoC'^) 



Fig. 3: The generating function F„ for rooted maps with root degree n is 
also that Zo^o(^) for Motzkin paths of length n with ordinate-dependent 
weights. 

Let us now come to the mobile interpretation of F^. From its definition of Sect. 1.2, 
it enumerates rooted maps with a root face of degree n. Taking as origin of the map the 
origin of the root edge, we may look at the clockwise sequence of distances from this 
origin of the successive vertices incident to the root face in the map. For convenience, 
in the planar representation, we choose as external face the root face itself so that the 
clockwise orientation around it corresponds in practice to the counterclockwise orien- 
tation around the rest of the map. Starting from the origin, this sequence of distances 
forms a Motzkin path of length n. The down-steps and the level-steps of the Moztkin 
path correspond respectively to the labeled vertices and the flagged vertices connected 
to the unlabeled vertex associated with the root face in the mobile. By removing this 
unlabeled vertex and all its incident edges, we obtain a sequence of rooted mobiles and 
half-mobiles (when removing an edge incident to a labeled vertex, we keep track of its 
position by distinguishing, say, the next edge incident to the same labeled vertex in 
clockwise direction). More precisely, on the Motzkin path, a down-step of the form 
(t, m) — 7- (t -(- 1, m — 1) corresponds to a rooted mobile with root label m satisfying (Rl) 
(possibly reduced to a single labeled vertex), while a level-step (f, m) — )■ (t + 1, m) corre- 
sponds to a half-mobile with root flag m also satisfying (Rl). Conversely, starting from 
a Motzkin path and a sequence of rooted and half- mobiles so associated with its down- 
and level-steps, a complete mobile is immediately obtained by connecting the mobiles 
to a new unlabeled vertex, and properties (P)-(R1)-(R2) are satisfled by construction. 
This clearly forms a bijection. Translated in the language of generating functions, it 
says that F^ is equal to the generating function for Motzkin paths of length n with 
ordinate-dependent weights {Rm)m>i and {Sm)m>o (see Fig. 3), i.e. F^ = ZQ^iji) with 
the notations of Sect. 2.1. Equation ( |1 . 1| ) follows from the general combinatorial theory 
of continued fractions. 
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To conclude this section, we mention that the continued fraction expansion may be 
directly interpreted at the level of maps, with no recourse to the coding by mobiles. It 
makes use of a particular decomposition of pointed rooted maps into slices, as explained 
in Appendix A. 

2.4- Related results 

Before proceeding to the next section, let us list a few other enumerative conse- 
quences of the coding of maps by mobiles. First, Rn and Sn satisfy recursive equations 
which translate decompositions of the corresponding mobiles. Indeed, let us consider 
a rooted (resp. half-) mobile with root label (resp. flag) n satisfying (Rl), not reduced 
to a single labeled vertex. We now decompose the mobile around the unlabeled vertex 
incident to the distinguished edge (resp. adjacent to the root flagged vertex). The three- 
step path associated with this unlabeled vertex has an arbitrary length k > 1 (thus leads 
to a weight gk), is positive but does not necessarily attain 0. The distinguished edge 
yields a distinguished down-step (resp. level-step) starting from ordinate n. Hence all 
the other steps form a three-step path of length k — 1 starting at ordinate n — 1 (resp. 
n) and ending at ordinate n. As in the above discussion of Fn, each down- or level-step 
starting at an ordinate m is associated with a rooted or half-mobile with root label 
or flag m satisfying (Rl) (except for the distinguished level-step when decomposing a 
half-mobile), and the decomposition is bijective. In the end, this translates into the 
relations 



Rn = l + '^gkRnZn-lAf^- 1) = 1 + ^ fi'fc^n.n-l (^^ " 1) 

fc=l 

oo 
J-n = / ^9kZn,n{k — 1) 



fc-1 fc-1 ^2.4) 



k=l 



where 1 stands for the mobile reduced to a single vertex labeled n, and the Z's are defined 
as in Sect. 2.1. In the second equation for Rn, we used the identity Zn,n-i{k — 1) = 
RnZn-i,n{k — 1) obtained by refiecting the paths with respect to a vertical axis. Note 
that the Z's depend implicitly on the sequences {Rni)m>i and {Sm)m>o hence relations 
( p.4| ) are of recursive nature. It is easily seen that they uniquely determine all R^s and 
Sn's as power series in the gkS. These equations were already derived in H and also in 



112011 using a different coding of maps (the bipartite case was first discussed in |jTT[] where 
an explicit solution was guessed). 

Let us now briefly discuss the simpler case of mobiles in which condition (Rl) is 
waived. By a simple shift of labels by — n, Rn (resp. 5"^) enumerates rooted mobiles 
(resp. half-mobiles) with root label (resp. flag) and with their labels all strictly larger 
than — n and their flags all larger than or equal to — n. Sending n — )■ oo amounts to 
waiving this lower bound. This implies that Rn and Sn converge for n — )■ oo in the 
sense of power series, i.e. each of their coefficients stabilizes. We denote by R and 
S their respective limits, which are nothing but generating functions for respectively 
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rooted mobiles with root label and half-mobiles with root flag 0. By the discussion 
at the end of the Sect. 2.2, we flnd that 2R + S'^ is the generating function for pointed 
rooted maps. From ( p.4| ) we deduce 



R = 1 + J29k\^P-iik - 1; R, S) 
k=i 

oo 

S = J29kP{k~l;R,S) 



(2.5) 



k=i 



with P_i and P defined as in Sect. 1. Note that these relations are tantamount to ( |1.6| ) 
by virtue of the identity P{k; R, S) = SP{k - 1; R, S) + 2VRP_i{k - 1; R, S). Note 
also that the k = 1 term in the first line of ( |2.5|) vanishes so that we may start the 
summation at k = 2 for R. 




external face of degree n 




(n,0) 



P{n;R,S) 



Fig. 4: The generating function F* for pointed rooted maps with root 
degree n is also that P(n; R, S) for three-step paths of length n starting 
and ending at ordinate 0, with ordinate-independent weights. 

We finally consider pointed rooted maps whose root face has a prescribed degree 
n, and denote by F* the associated generating function (where the root face does not 
receive a weight Qn)- F* only differs from F^ by the extra marking of a vertex or, said 
otherwise, we no longer impose that the origin of the map be the starting point of the 
root edge. Denoting by d the distance from the origin of the map to the origin of the 
root edge, we now label each vertex by its distance from the origin of the map minus d. 
Applying the mobile construction rules of Fig. 1, we obtain a mobile with a distinguished 
unlabeled vertex corresponding to the root face. Its associated three-step path codes 
the distances (minus d) from the origin to the map vertices incident to the root face. 
The origin of the root edge distinguishes a step starting at ordinate on the path, 
therefore the path can be seen as a three-step path of length n which starts and ends at 
ordinate 0, but which is not necessarily positive. Decomposing around the distinguished 
unlabeled vertex, each down- or level-step of the path yields a rooted or half-mobile, 
whose root label or fiag we may set to by a suitable shift. This decomposition is 
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bijective, and translates into the relation (see Fig. 4) 

F:^P{n-R,S). (2.6) 

This relation is analogous to the relation F^ = ZQ^Q{n) established above, yet it is 
considerably simpler because it involves three-step paths with weights independent from 
the ordinate. In this sense, it is of the same nature as formula ( p..4| ) , which we shall now 
prove. 

3. An explicit expression for the moments 



We now turn to the derivation of formula ( |1.4| ) for the moments Fn as defined in 
Sect. 1. In contrast with the J-fraction expansion ( |1 . 1| ) , which follows from a general 
scheme where planar maps and mobiles turn out to fit, formula (|l]^) seems deeply 
related to the map structure. One of our purposes is to show that it naturally results 
from several approaches. First, we discuss a direct yet non-constructive check based on 
the comparison with formula ( p.6|) (Sect. 3.1). Then, we explain how formula ( |1.4| ) arises 
from the solution of Tutte's equation, which is the original approach to enumerating 
maps (Sect. 3.2). Finally we present a bijective proof involving distance-dependent 
generating functions and mobiles (Sect. 3.3). Beyond our main purpose, this third 
approach yields a one-parameter family of expressions for F^, the so-called "conserved 
quantities" , interesting on their own. 

3.1. Proof via pointed rooted maps 



It has been noted on several occasions |]21[] that considering maps that are both 



pointed and rooted brings some simplifications. This is illustrated by formula (|2.6| ) in 
our context. A natural idea is to deduce formula ( |1.4| ) from it. Because F* only differs 
from Fn by the extra marking of a vertex, there is a direct relation between these two 
quantities. 

To make this relation explicit, it is most convenient to introduce an extra weight 
u per vertex. We denote by F*{u) and Fn{u) the correspondingly modified generating 
functions. (Note that u is actually a redundant parameter, as by a simple counting 
argument and Euler's relation we have Fn(u) = u'^/'^+^ Fn(l)l ,„ ,,k/2-i and the same 
for F*.) On the one hand, relation ( p.6|) immediately generalizes as 



F^{u)=uP{n;R{u),S{u)). (3.1) 

where R{u) and S{u) enumerate mobiles and half-mobiles with an additional weight u 
per labeled vertex, and the extra factor u accounts for the origin. On the other hand, 
the extra marking of a vertex yields 

F:{u)=uj^Fn{u). (3.2) 
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Therefore Fniu) is a primitive of P{n; R{u), S{u)) with respect to the variable u. It is 
shown in Appendix B that it has the exphcit form 

oo 

F,(w) = Y, Aq{u)P+{n + q- R{u),S{u)) (3.3) 

with Aq{u) defined as in ( |L5|) with R,S replaced by R{u),S{u). By specializing to 
u = 1 we recover the expression ( |1.4| ) for F^ = -Fn(l)- 

3.2. Tutte's equation and its solution 

The generating function F^ for rooted planar maps with root degree n is of ut- 
most importance in Tutte's original approach 0. Indeed, the whole family {Fn)n>i is 
uniquely determined as a power series in the variables {gk)k>o by the equation 

n — 2 oo 

Fn = 2_^ FiFn-2-i + 2_^9kFn+k-2 (3.4) 

i=0 k=l 

valid for all n > 1, with the convention Fq = 1. This equation directly expresses that, 
removing the root edge of a planar map with root degree n, two situations may arise: 
either the map is split into two connected components, which can be seen as two rooted 
planar maps whose root degrees add up to n — 2 (possibly one of these maps is reduced 
to a single vertex and has root degree 0), or the map is not split, then its root degree 
is increased hy k — 2 where k is the degree of the face formerly on the left of the root 
edge. Passing to F{z) we get Tutte's equation 

oo / fc-2 

F{z) = 1 + z^F{zY + Y. ^^^'"' ^(^) - E ^'^^- I (3.5) 

k=l \ j=0 

also known as loop equation in the context of matrix models. 



The solution of Tutte's equation ( p.5|) was obtained by Bender and Canfield |J22 
(with slight restrictions on the Qk), see also [^^ for a more general formulation and |^ 



for the matrix model counterpart. These authors were ultimately interested in the series 
F2 which, upon squeezing the bivalent root face, yields the "true" generating function 
for rooted planar maps. Our purpose here is to show that the general formula ( p..4| ) for 
Fn follows straightforwardly from their discussion, which we now recall briefly. 

We first observe that Eq. ( p.5| ) is quadratic in F{z) (viewing the terms Xl^^o ^^ -^j 
as constants), hence F[z) is readily given by 



nz) = ^ 1 - E 3kz^~' - ^^) (3-6) 



fe=i 
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where A{z) = (l - J2T=i 9kZ^ '') ^ 42;^ (l - J2T=i 9k Yj)=1 ^^ ^^^ Fj\ is the discrim- 
inant. As shown in [ p2l - p3| , ^J /S.{z) has a factorization of the form 

oo 

^/X{z)=T{z-^)^/^, r(z-i) = ^7gz-^ t^{z) = I + ti^z + K2Z^ (3.7) 

where the coefficients 7^, ki, K2 are power series in the gt to be determined. (In 
the actual proof, one must for a while assume that degrees are bounded, i.e. that gk 
vanishes for all k larger than some fixed integer, which from (^) implies that vA( 



does not contain arbitrarily large negative powers of z. In this situation, the existence 
of a factorization follows from Brown's theorem [^^, and T{z~^) is a polynomial in 
z~^ . The bound on degrees can be lifted eventually and the product T{z~^)^/k{z) has 
a well-defined expansion in z.) 

Let us now explain how to compute practically k(z) and T{z~^). Following the 
approach of 0, we introduce the power series R and S defined by 

5 = -^, i? = ^i^i^, i.e. K{z) = {l-Szf-ARz''. (3.8) 
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This choice is particularly convenient because the series expansions of y Hi[z) and 



1/ \J k{z) are related to the generating functions P(n; R, S) and P^{n; R, S) for three- 
step paths defined in Sect. 1. Namely, it is elementary to check that 

00 ^00 
\^<iy=l-^2-2i?22^P+(n;i?,^)z", = ^P(n;P,g)z". (3.9) 

n=0 V^(^) n=0 

Now, substituting (|3.7|) in (|3.6|) , extracting the coefficient of z"^ for n > and using the 
first relation ( |3.9p , we obtain 

00 
F^ = RY,l<iP^in + q;R,S) (3.10) 

q=0 



which coincides with (|1.4| ) up to the identification Ag = ■jqR. It remains to determine 



the unknowns 7^, i?, S. For this it is convenient to rewrite ( |3.6| ) as 



00 



z-^-E dkz-'' - 'iF{z) 

and recall that F{z) only contains nonnegative powers of z, with Fq=1. Using then the 
second relation (|3.9[ ), extracting the coefficient of z~'^~'^ for g > leads directly to ( |1.5| ) 



while extracting the coefficients of z ^ and z^ yields nothing but (|1.6| ) . This establishes 
the general formula for F^ announced in Sect. 1. 
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d{Vo,v)^d 

d{Vo,v) = J (i^o, boundary) 




(n,d) 



z;/«) 



Fig. 5: The generating function F^-d foi' pointed rooted maps with root 
degree n, where the origin v of the root edge is one of the vertices closest 
to the origin vq of the map among those incident to the root face, with 
a distance d{vo,v) less than d, is also that Z^^{n) for three-step paths of 
length n starting at, ending at, and restricted to stay above ordinate d, 
with ordinate-dependent weights. 



3.3. Bijective derivation 

Let us finally present a last derivation of the expression (|l]^) for the moments F^ . 
It is based on a bijective decomposition of pointed rooted maps keeping track of some 
distances between the origin and vertices incident to the root face. Passing to mobiles, 
this gives rise to a one-parameter family of expressions for F^ (the so-called conserved 
quantities) , interesting on their own, and encompassing ( p..4|) -( pr5D as a limit case. 

We first introduce the generating function fn-i {i > 0) for pointed rooted maps 
whose root face has degree n and whose origin is both at distance i from the origin 
of the root edge and at distance i from the boundary of the root face, defined as the 
set of vertices incident to this root face. In other words, we demand that the origin of 
the root edge be one of the vertices closest to the origin of the map among all vertices 
incident to the root face. For i = 0, the origin of the map is forced to be the origin of 
the root edge, so that we have simply /^;o = F-n- We also introduce Fn-^a = J2i=o fn;i 
enumerating maps where the above distance i is less than or equal to some fixed value 
d. Clearly, we have for any d > the equality 



n:d 



E/" 



(3.12) 



The two terms in the r.h.s may be evaluated in the mobile language as follows: to 
evaluate Fn-^d, we look at the sequence of distances from the origin of the vertices 
around the root face. Starting from the origin of the root edge, we have a three-step 
path from (0, i) to [n, i) for some i between and d, with all its ordinates larger than or 
equal to i. This path defines a sequence of labels and fiags for the labeled and fiagged 
vertices around the unlabeled vertex associated with the root face and the full mobile 
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is obtained by attaching a rooted mobile to each of these labeled vertices and a half- 
mobile to each of these flagged vertices. The entire mobile must satisfy (Rl) and (R2) 
but, by the same argument as in Sect. 2.3, we may transform it into a mobile satisfying 
(Rl) only by shifting all labels by d — i. The attached rooted mobiles (respectively 
half-mobiles) are now counted by Rm (respectively Sm) with indices corresponding to 
the down- (respectively level-) steps of a shifted three-step path from (0, d) to (n, d) 
which never dips below d. To summarize, Fn-d is the generating function of three- 
step paths from (0, d) to (n, d), staying above d, and with a weight Rm per down-step 
(t, m) — )■ {t + l,m — 1) and Sm per level-step (t, m) — )■ (t -(- 1, m) (see Fig. 5 for an 
illustration). With the notations of Sect. 2.1, our first result may be stated as: 



n\d 



ztM 



(3.13) 




{n,i) 




(n,i+j) 



(n+k-l,i-l) 



Fig. 6: A schematic explanation for the bijection between maps in f^-i for 
i > 1 and maps in Ai{n;i, k,j) with arbitrary /c > 3 and j > 1 (see text). 



As for the subtracted term in ( p.l2| ) , it concerns maps where the origin of the root 
edge, which we from now on denote by v, is at distance i > 1 from the origin of the 
map denoted by vq. Then there exist vertices adjacent to v and at distance i — 1 from 
vq, hence edges of the type {i,i — 1) incident to v. Let us pick the leftmost such edge 
(ordering them from the root edge, see Fig. 6) and denote by v' its endpoint. Note that 
v' cannot lie on the boundary of the root face as, otherwise, the distance from vq to the 
boundary would be strictly less than i. Therefore, the face on the left of the leftmost 
edge from v to v' cannot be the root face but is instead a regular face of degree k > 3 
weighted by gk (the degree cannot be 2 since we picked the leftmost edge from t; to f'. 



19 



nor 1 since it is of type (z — 1, z)). We may now merge this face with the root face into 
a new external face with larger degree n + k hy simply splitting the vertex v into two 
vertices v and v as shown in Fig. 6. Looking again at the sequence of distances along 
this larger boundary starting from the new vertex v, the first n + 1 vertices correspond 
to vertices of the former boundary and their distance to vq is necessarily larger than 
or equal to what it was before since the splitting did not create new adjacences but 
suppressed some. Moreover, the new vertex v is still at distance i from vq but v is 
necessarily at a distance strictly larger than i. This is again because we chose the 
leftmost occurrence of an edge (i, i — I) adjacent to v. The sequence of distances, when 
going from t; to -y around the new external face, now forms a three-step path from (0, i) 
to (n, z + j) for some j > 1, that never dips below i. The following distances (going 
from V to v') then form a three-step path from (n, i + j) to (n + k — l,i — 1) with no 
imposed lower bound (apart from having non-negative ordinates), as displayed in Fig. 6. 
To summarize, we end up with a pointed rooted map, whose root face has degree n + k 
and whose sequence of vertex distances around it defines a three-step path made of 
three parts: a path from (0, z) to (n, z + j) for some j > 1, that never dips below z, a 
three-step path from (n, i + j) to (n + /c — 1, z — 1) that never dips below and a final 
up-step (n + /c — 1, i — 1) — )■ {n + k,i). We call Ai{n; i, /c, j) the family of such maps. 



If conversely, we start from a map in Ai{n]i, k,j) for fixed n > and i > 1 and 
arbitrary A; > 3 and j > 1, gluing the first vertex v of the boundary to the n-th one 
V into a single vertex v will split the external face into a face of degree k and a new 
external face of degree n. The new vertex v will be at distance i from the origin and all 
the vertices on the new boundary will be at distance larger than or equal to i. Indeed, 
the only possible way to reduce the distance from the origin after gluing is to use a 
path which goes through v but such a path has a length larger than or equal to z so 
the distance cannot be reduced to less than z. Finally, all the vertices adjacent to v 
were at a distance larger than or equal to z + j — 1 > z hence their distance after gluing 
remains > i so that the leftmost edge {i,i — 1) incident to v after gluing is precisely 
the edge from v to the vertex v' formerly preceding v along the external face. In other 
words, the splitting and gluing transformations are inverse of one-another and provide 
a bijection between maps enumerated by fn;i for some fixed n > and i > 1 and maps 
in Ai{n;i,k,j) with arbitrary k > 3 and j > 1. These latter maps may themselves 
be transformed into mobiles satisfying (Rl) and (R2). As before, the set of all values 
i = l,...,(i is well captured in a single family of mobiles satisfying (Rl) only upon 
shifting the labels and flags hy d — i. In these mobiles, the three-step path ensuring 
(P) around the unlabeled vertex associated with the external face is now a path from 
(0, d) to (n -|- /c, d) made of three parts: a first part of length n from (0, d) to (n, d + j) 
for some J > 1 that never dips below (i, enumerated by Z'j"^^(n), a second part of 
length k — 1 from (n, d + j) to {n + k — l,d — 1) that never dips below 0, enumerated 
by Zd+j,d-i{k — 1), and a final up step from {n + k — l,d— 1) to {n + k,d), with weight 
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1. This is summarized in the equahty: 

d 

J2fr^-^^ = Y.9i^Y.^h+M)Za+j,a-i{k - 1) . (3.14) 

i=l fc>3 j>l 

From (|3.13| ) and ( |3.14| ), we deduce eventuaUy the expression for F^: 

Fn = ZX,{n) -Y,9kY.ZX,^M)Za+j,a-i{k - 1) (3.15) 

A;>3 j>l 

valid for any d > 0. An important outcome of this result is that the quantity on the r.h.s 
of Eq. (|3.15|) is a conserved quantity, i.e. it is independent of d. This conservation prop- 
erty was aheady proved in [53 in the restricted case of bipartite maps via completely 



different and much more involved arguments, with no direct combinatorial interpreta- 
tion. Taking d = 0, we have in particular the equality F^ — Z(^Q(n) which corresponds 
precisely to the interpretation of F^ which was discussed in Sect. 2.2 and which led us 
to the continued fraction formula ( |1 . 1|) for F{z). Another interesting limiting case is 
when (i — )■ oo in which case the r.h.s in ( |3.15|) may be expressed in terms of R and 5" 
only. Before we take this limit, let us derive two slightly different, although essentially 
similar, expressions for F^. They can be obtained by recalling the second recursion 
relation (|2.4|) for Rd, which be may used to rewrite ( |3.15| ) in the slightly different form 



Fn = Z+,{n)Ra -J^^'^Yl ^i+,(^)^d+,,d-i(/c - 1) (3.16) 

k>2 j>0 

where the sums now include the contribution of j = and k = 2 (note that the term 
/c = 1 in the first line of (|2.4|) is zero). The enumerated three-step paths are now paths 
from (0, d) to {n + k — l,d—l) whose first part of length n has ordinates larger than or 
equal to d, with the ordinate d now allowed for the abscissa n. Denoting hy n + q + 1 
the abscissa of the first occurrence of the ordinate d — 1, with q between and /c — 2, 
we may write alternatively 

Fr, = Z+^in)Rd -Y.9kY. ^Xdi^ + <i) Rd Zd-i,d-iik -2-q) (3.17) 

fc>2 q>0 

valid for all d > 0. In particular, by letting d — )■ cxo and recalling from Sect. 2.1 that 
Zd,d{n) -^ P{n; R, S) and ZJ~^(n) — )■ P~^{n; R, S), we arrive at the simplest formula 

Fr, = P+{n; R, S)R - ^ £?fc ^ P+(n + g; R, S)RP{k-2- q; R, S) (3.18) 

k>2 q>0 

which is precisely the desired expression ( p..4| )-( |T75D . 

A slightly more direct derivation of ( p.l8|) consists in using mobiles without the 
conditions (Rl) and (R2) as done at the very end of Sect. 2.2. More precisely, we may 
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d(Vo,v) =<i (1^0, boundary) 
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P(n;R,S) 



Fig. 7: The generating function Xli^o ■^"^;* ^°^ pointed rooted maps with 
root degree n, where the origin v of the root edge is one of the vertices clos- 
est to the origin vq of the map among those incident to the root face, is also 
that P~^{n; R, S) for Motzkin paths of length n with ordinate-independent 
weights. 

in the first place label the vertices of the maps counted by fn;i by their distance to 
the origin minus i and consider values of i no longer bounded by d. We end up with 
mobiles where the conditions (Rl) and (R2) are waived and whose sequence of labels 
and flags around the unlabeled vertex associated with the root face match the down- 
and level-steps of a three-step path from (0, 0) to (n, 0) that never dips below 0. In the 
absence of the conditions (Rl) and (R2), each down- (respectively level-) step is now 
weighted by R (respectively 5") so that we arrive directly at 



Y,fn■^ = P+{n■,R,S). 



(3.19) 



This formula is illustrated in Fig. 7. If we now consider the case of arbitrary i > 1 and 
use the above bijection of Fig. 6 with configurations of rooted maps in Ai{n]i, k,j) for 
arbitrary k > 3 and j > 1, we arrive by again shifting the labels in these maps by i to 
mobiles counted directly by 



/ ^ ■ln;i 



(3.20) 



J2 9k J] v^P+(n; R, 5)P-,-i(A; - 1; i?, S) 

i=l A;>3 j>l 

where, for an arbitrary integer to, PmiQ', R^ S) denotes as before the generating function 
for arbitrary three-step paths from (0, 0) to (g, to) with a weight yR per up- or down- 
step and S per level-step and, for to > 0, P^{n; R, S) restricts the enumeration to paths 
whose ordinates all stay above 0. Using the identity ( p.5|) for R we arrive directly at 

k>2 j>0 

I- R, S)R - ^ £?fc ^ P+(n + q; R, S) RP{k 

k>2 q>0 

upon denoting, in the r.h.s of the second line, hj n + q + 1 the abscissa of the first 
occurrence of the ordinate —1 in the paths enumerated in the r.h.s of the first line. 



Fr, = P+in;R,S)R 



>+, 



n; 



RP^in; R, S)P-j-i{k - 1; R, S) 

2-q-R,S) 



(3.21) 
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4. Study of the Hankel determinants 



Having establistied tlie formula (|1.4|) for F^, we now study its consequences on the 
Hankel determinants Hn = det Fi+j and minors Hn = det Fi+,+5 . We first 

0<i,j<n ■" 0<i,j<n ■' ''' 

establish formulas ( p. .71) and ( |1.8| ) via a natural path decomposition which translates 
into matrix identities (Sect. 4.1). We then identify the resulting determinants with 
symplectic Schur functions (Sect. 4.2). 



4^.1. Path decomposition 



P^_l{q;R,S)-P^^l^2i^;R,S) 



P^(i;R,S) 



•&. 



P (i+j+q,R,S) 



P, (j;R,s) 




Fig. 8: A schematic picture of the path decomposition leading to (4.1). 
We first establish the formulas (|1 . 7|) and ( |1.8| ). We start with the identity 



P+{i+j + q;R,S) = J2J2^ki^-^R^S) {Pk-e{q;R,S)-Pk+i+2{q;R,S))P+{j;R,S) 

k=o e=o 

(4.1) 

obtained by decomposing the paths of length i + j + q in P+(z + j + q; -R, 5") into a 

first path of length i ending at the ordinate k, a second path of length q ending at the 

ordinate i and a third path of length j back to ordinate (see Fig. 8 for an illustration) . 

The medial part must have non-negative ordinates only, hence we must subtract from 

Pk-e{q; -R, S) the configurations which dip below 0. By a classical reflection argument, 
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those configurations are in one-to-one correspondence with paths from ordinate —k — 2 
to ordinate £, enumerated by Pk+£+2{Q', -R, S). 

Using now formula ( |1.4| ), we deduce the matrix identity: 

H = T^ ■ B ■ T (4.2) 

where H, T and B are the semi-infinite matrices 

H = (-Pi+j)i,j>o, T = (P^ {j;R,S))ij>o, B = {Bi^j — Bi+j+2)i,j>0i (4.3) 

with Bi defined as in j ^TT^ ) and T""" denoting the transpose of T. Note that T is upper 
triangular, hence the matrix product ( ^.2| ) involves only a finite sum for each matrix 
element of H. 

For a semi-infinite matrix M = {Mij)ij>o, let us denote by M„ its submatrix 
obtained by keeping only the row and column indices between and n. Note that M„ 
has size {n + 1) x (n + 1). We also denote by M^ , M^ and M^'"' the submatrices 
obtained by further removing respectively the row with index i, the column with index 
j and both row i and column j. 

We have, by definition, H^ = detH^. Because T is upper triangular, the identity 
( ^?2|) restricts safely to 

H^ = T^ ■ B„ ■ T^ (4.4) 

which yields 

i7^ = (detT^)^detB^. (4.5) 

Noting that T^ = P^{i) = -R*/^, the triangular determinant is readily computed, estab- 
lishing formula ( [1.7|) . 

As for Hn = detH^";|_i '" , we have the restriction 

TT{n+l;n) _ rpT r,{n+l;) rp{;n) / , ^x 

■"■n+1 — -"-n ' ^n+1 ' ^ n+1 l^-"J 

and the Cauchy-Binet formula yields 

n+l 

H^ = det T„ Y. det Bi",Y^^^ det T^:;"^ (4.7) 

i=0 

Note that T^^^'J^i is still triangular and, unless i = n or n + 1, it has at least one 
zero on the diagonal. Therefore only those two terms are non-zero in ( |4.7| ) and, as 
Tn,n+i = Pni''^ + 1) = (^ + 1)R^^'^S, the triangular determinants are again readily 
computed. Finally we have det bJ;^+/'''+^^ = detB„ = if^/i^'^^'^+i)/^, and detBj;^^'''^ 
is nothing but the determinant in formula (|1.8|) , which is now established. 



To conclude this subsection, let us mention that using the decomposition ( [4.2| ) of 
the Hankel matrix H for determinant evaluation amounts to doing elementary row and 
column manipulations, as stated in Sect. 1. 
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4-2. Determinants and minors of the matrix B identified as Schur functions 

We now turn to the study of the quantities det B„ and det B^ '" encountered 
above, which are minors of the semi-infinite matrix B. This matrix is remarkably the 
difference between the symmetric Toephtz matrix (i?i_j)i j>o and the Hankel matrix 
(i?i_i_j_i_2)i,j>05 strongly suggesting possible determinantal identities. The Laurent series 
Si^-oo -^j^* i^ known as the symbol of the Toeplitz matrix. 

To proceed further with our analysis, we shall assume that the map degrees are 
bounded, namely we fix a positive integer p and consider maps whose faces have degree 
at most p + 2. As already mentioned above, this amounts to setting gk = for k > p + 2, 
so that we only deal with power series in the finite set of variables (71, ... , gp+2- Then Bi 
vanishes for \i\ > p (and Bp 7^ 0) so that B is a band matrix. The symbol is a symmetric 
Laurent polynomial of degree p in x, which thus admits 2p roots of the form x = 
(xi, l/xi, X2, 1/^2, . . . , Xp, l/xp). These are precisely the solutions of Eq. (|1.9|) (which 
by analogy with the "perturbative" approach of |]11[ we call the characteristic equation) . 
Up to the overall factor Bp = —gp^2R, we may view x as a parametrization of 
the S's, namely 

B, = {-ir+^Bpep+,i^) (4.8) 

where e^ (x) is the ith elementary symmetric function of x, defined for instance via 



P . X 00 

l[{l + x,t)il + -) = J2 e^w^^- 

1=1 ^ ' i = — oo 



(4.9) 



Note that ep_j(x) = ep-i-j(x) and ei(x) = for z > 2p or z < 0. It follows that detB„ 

and det B^ '" are symmetric functions of x, which we may hopefully identify. 

Let us first quickly describe a naive strategy, based on finding the zeroes of the 
symmetric functions at hand. We observe that, for x a root of the symbol, the semi- 
infinite vector ((x*"*"^ — x~'^~^)/{x — x~^)) .^ is in the kernel of B. Its restriction to 

the indices {0, . . . , n} is however generally not in the kernel of B^ nor Bn '^ ■ We 
may attempt to find a vector in the kernel of B^ (resp. B^ '^ ) by taking a linear 
combination over the p independent roots of the symbol. Doing so, we find that we have 
to satisfy p "boundary conditions" which are linear equations for the p coefficients of 
the linear combination. A non-trivial solution exists if the associated p x p determinant 
vanishes, hence this determinant (which is again a symmetric function of x) divides 
detB„ (resp. detBn '^ )■ By comparing degrees, we find that they only differ by a 
constant. In the end, we have expressions for det B„ and det B„ '" in terms of p x p 
determinants. It turns out that these are nothing but the determinants present in the 
Weyl character formula for the symplectic group Sp2p, written down below. What we 
have just encountered are classical formulas in representation theory. 

Let us now be more educated and directly identify det B„ and det B^T '"^ as 
instances of the general "symplectic e-formula" |1T7|,|T8| 



sp2p(A,x) = det (ex'-j+ii^) - eA'_j-j(x) . (4.10) 
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Here sp2p(A, x) stands for the symplectic Schur function associated with the partition A 
(having at most p parts), A' denoting its conjugate partition (having at most m parts). 
Comparing with ( |4.8| ) and recalhng that ep_|_g(x) = ep_g(x) for aU q, we obtain the 
general formula 

det (S.-,-^, - S.+,+M,+2) = (-l)l^li?;+isp2p(A,x) (4.11) 

0<.i,j<.n 

where fXj = p— A^_|_j^, m — n + 1 and |A| denotes the sum of A (equal to that of A'). It is 
now immediate to identify det B^ with the case /ij — 0, hence with the "rectangular" 
partition Xp^^+i = p'^+^ made of n + 1 parts of size p, conjugate to the partition 

Ap,n+i = {n + 1)P with p parts of size n + 1. Similarly detB^ '"^ corresponds to the 
case Hj = Sj^n, hence to the "nearly-rectangular" partition X'p^n+i = p^{p — l) made of n 
parts of size p and one part of size p — l, conjugate to the partition Ap^^+i = (n + l)P~^n 
with p— 1 parts of size n + 1 and one part of size n. We end up with the compact formula 
for the Hankel determinants 

H^ = {-ir^^+'^ B;+' r""^^ sp2p(Ap,^+i, x) (4.12) 

from which, using ( |1.2|) , we deduce 

„ O^P2p(Vn+l,x)sp2p(Ap,^_l,x) 

H,, = K — — — ^ . (4.1d) 

Similarly, for the Hankel minors we have 

H^ - {n + l)SH^ = (-l)p(-+i)+is;+ii?^^^^^ sp2p(Ap,„+i,x) (4.14) 

from which, using ( |1.3|) , we deduce 

(4.15) 



g -g /^/^ SP2p(Ap,^+l,x) SP2p(Ap,n,x) 

SP2p(-^P,n+l7x) SP2p(Ap^n,x) 



Beside the above e-formula, other expressions are known for the symplectic Schur 

function sp2p, namely the "/i- formula" and the Weyl character formula. Both involve 

determinants of size p, independently of the variable n in Hn, Hn, Rn and Sn- On the 

one hand, the /i- formula involves the complete symmetric function /ii(x), defined for 

instance via 

P -I -I °° 

^^37^T3x = E^^w^^ (4-16) 

i=l * Xi i=o 



and reads [17-18 



SP2p(-^.x) = det ( /lA,._j + l(x) ; /lA,_j+z(x) + /lA,_j-i+2(x) ) (4.17) 
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where [aj [cLij) denotes the matrix with elements aj in the first row and elements aij 
in the rows i > 1. It immediately allows to rewrite the r.h.s of ( 4.12 ) and ( [4.14| ) as 



determinants of size p. On the other hand, the Weyl character formula reads [^ 

det {x^''+'~'+' -x~^''+'~'~') 

sP2p(A, x) = '-^■^-" ^ ^ ^ ^. — (4.18) 

det [x^ — Xj J 

and involves the ratio of two determinants of size p. However, the denominators cancel 
in the expressions ( [4.1tj| ) for R^ and ( |4.15D for 5"^, which establishes the announced nice 
formulas ([TTlD and ( pH^ ). 



5. The special case of maps with even face degrees 

As mentioned in Sect. 1, some simplifications occur when we restrict our analysis 
to bipartite maps. The distance from a given origin in the map changes parity between 
adjacent vertices and consequently, there are no edges of type (n, n), leading to mobiles 
without flagged vertices. At the level of generating function, this implies that 5" and Sn 
for all n vanish so that the continued fraction ( |1 . 1| ) becomes of the Stieltjes type ( |1.13D , 



consistent with the property that F^ = for n odd as the root face must have even 
degree in a bipartite map. The Hankel determinants factorize as 

H2n = hl^^hi% H,^+, = h(^^h^^^ (5.1) 

where 

/^f = , det i^2.+2„ hl^^ = det F2.+2,+2 ■ (5.2) 

Using the same path decomposition as in Sect. 4.1, we may write 



U<.i,j<-n 

hW = Rin+lf det (S,_, - S,+,+2) 
U<.i,j<n 



(5.3) 



where we introduce the notation 

B, = B2^ (5.4) 

for all integers i (note that -622+1 = and -B_i = Bi). Let us now turn to the case 
of maps with maximal degree 2p + 2, for some p > 1, so that Bi vanishes for \i\ > p. 
Writing the characteristic equation as 

p 
J]S,y* = 0, y = x^, (5.5) 
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whose 2p roots are gathered in the 2p-uple y = (yi, l/yi, . . . ,yp, l/Vp), we now have the 
identification 

B, = {-ir+^Bpep+j{y) (5.6) 

in terms of the elementary symmetric function ej(y). This yields 

h(0) = (_i)P(^+i) (S,)-+ii2-(-+i) det (e,_,+,(y) + e,_,_,+i(y)) 

= (_l)p(n+l) (Sp)n+l^n(n+l) q^^+i (Ap,^+i, y) 

;,« ^ (_i)P(n+i) (s^)n+i^(n+i)^ ^^^ (ep-,+.(y) - e,_,_, (y)) ^^'^^ 

= (_1)P(-+1) (S,)-+ii?("+i)' sp2p(Ap,.+i,y) 

where, as before, Xp^n+i = (n + 1)^ is the partition with p parts of size n + 1. Here we 
used again ([4 .101) to identify the second determinant with a symplectic Schur function 



of y, while the first determinant is now recognized as an instance of the general formula 



02p+i(A,y) = det eA'-j+i(y) ;eA'-j+z(y) + ca' j-^+2(y 

l<-i,j<.m \ ■' ■' ■' 

= , det fex' j+i{y) + e^' j-i+i{y 



where 02p+i(A,y) stands for the odd- orthogonal Schur function associated with the 
partition A (having at most p parts). A' denoting its conjugate partition (having at 
most m parts), and where y denotes the (2p + l)-uple (yi, l/j/i, . . . ,yp, l/yp, 1) with an 
additional 1 term. The first identity may be found in fl^ and it implies the second 
one by elementary manipulations on the determinant. Again, the odd-orthogonal Schur 
function admits a simple Weyl formula, namely: 

02p+i(A,y) = ~''^~^ — T — ^T ■ (5.9) 

, det {y^ ^ -y^ ' 



2 



1^ ' (~i-m/H Ai^ -^ 



Alternatively, the factorization (|5.1|) of H^ , with hn and hn expressed directly as in 
( p.7| ) , may be obtained right away from the general formula ( [4.12| ) (with p replaced by 
2p) for Hn upon using the identity 

sp4p(A2p,n+i; X U -x) = (-l)P('"+i) sp2p(Ap^ ^n+i J , y) 02p+i(Ap_ ^^j , y) (5.10) 

where x U — x = (xi, 1/xi, . . . , Xp, 1/xp, — xi, — 1/xi, ... — Xp, — 1/xp) and where y = 
(yi, 1/yi, . . . , l/p, 1/yp) with yi = (x^)^. This identity is easily proved via elementary 
determinant manipulations in the Weyl formula for each of these Schur functions. 
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The final "nice" formula for R„ becomes 



Rn — R 



det 



.f+J+- 



-2-J- 



^ det [y^ 



t+j-i 



Vi 



-t-J+i 



det [yt'-y~^ ' 



det (yp-^-yr^"'^^) 

l<i,j<p v / 



(5.11) 



which matches formula (5.8) of Ref. [^, given there without a complete proof. To make 
the identification complete, we have to rewrite the characteristic equation (|1.10| ) as 



P+i fc-i 
l = 5Z^2fc5^P(2A;-2-2Q;i?,0) [VRx + 

k=l q=0 

p+1 k-1 



R 



2q 



fc=l g=0 ^ 

P+1 fc-l ^ 

E *'«'-' E ( 

k=\ m=0 ^ 



2k-2-2q\ 1 



2g 



2A;-1 
k — m — 



1. S 



j = -m 



(5.12) 



with y = x^. This is precisely the form found in |]Tl|]. The third step in ( |5.12| ) follows 
from the identity 



fc-i 

E 

'7=1 J I 



2k-2q-2\f 2q 
k-q-1 )\q+j 



k-i 

E 

m=\j\ 



2k -1 
k — m — 1 



(5.13) 



which itself follows from the more obvious identity 



fc-i 

E 



9=1^1 



2k-2q-2\ r/ 2q 
k-q-1 ) [\q+\j\ 



2q 

9 + lil + i 



2A;-1 

k-\j\-l 



(5.14) 



obtained for instance by enumerating paths with ±1 steps of length 2k — 1 starting 
at height and ending at height 2|j| + 1 in two manners: either directly (r.h.s) or by 
decomposing at the last passage at height (l.h.s). 



6. Triangulations and quadr angulations 

We now turn to the particularly simple cases of triangulations and quadrangula- 
tions, where the expressions ( p. . 1 1| ) , (|1.12| ) (for triangulations) and ( |5.11j ) (for quadran- 
gulations) involve only determinants of size 1. They admit an elementary combinatorial 
interpretation in terms of one- dimensional hard dimers, as we shall see. 
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-1-2-1 12 3 



Fig. 9: Interpretation of formula (6.3) for -Fi+j in tlie case of triangula- 
tions. The weights per step are different in the central strip as shown. 




-n 



-2-1 12 3 



n+i 



Fig. 10: Correspondence between configurations of non-intersecting paths 
on the graph of Fig. 9 and hard dimers on a segment of length n. 

6.1. Triangulations 

The case of triangulations, characterized by gu = g3Sk,3 is particularly simple. In 
this case, the only non- vanishing Aq coefficients are Aq and Ai, given from ( |1.5| ) by 

Ao = R-gsRS, A^ = -g^R (6.1) 

with R and S given implicitly from (|1.6|) by 



S = g^{S^ + 2R), R=i^^{^S^^%g^RS)-^ = l + 2gr,RS 
From (|1 . 7] ) , the non- vanishing Bi coefficients are 

Bo = AQ + AiS = l-2g^RS = l 
Bi = B_i=Aiy/R = -gs R^^^ 



(6.2) 



(6.3) 
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while the characteristic equation ( pTOl ) reads 

Formulas (|1.11|) and (|1.12P for the two-point functions reduce to 



(6.4) 



_ (x"+2 - a;-'^-2)(x" - x-^) 



5„ = >5 - v^ 



™n+l ™ — Ti— 1 



^71+2 ^ — 71—2 ^n+1 ^— n — 1 



(6.5) 



To recover the expressions given in [^, we set 

y = x^ , (6.6) 



which is solution of 



and obtain eventually 



1 

y ' " giR' 



2^ + 7 + 2 = -^—, (6.7) 



Rr, — R 



(1 -;/")(! -y"+2) 

(1-|/-+1)2 



5. = ^-,3i?V (1-^^(1-^^) 



(6.^ 



(l-y"+i)(l- 2/^^+2) ■ 



Interestingly, these formulas admit a combinatorial interpretation in terms of one- 
dimensional hard-dimers. Indeed, for triangulations, formula ( p..4| ) for the generating 
function Fn reduces to 

Fr, = Ao P+{n; R, S) + Ai P+{n + 1; R, S) (6.9) 

so that -Fi+j may be interpreted as enumerating three-step paths with non-negative 
ordinates from, say, (— ^,0) to {j + 1,0) with a weight a/R per up- or down-step and 
5" per level-step, except for steps in a central strip of width 1 (i.e. the strip between 
abscissas and 1) which receive instead weights incorporating the Ag factors, namely 
Ai X y/R — —Qs i?^/2 fQj. ii^Q up_ and down-steps, and AqxI + AixS = R — 2g^ RS = 1 
for the level-steps (see Fig. 9). Note that these weights are nothing but Bi and Bq 
respectively, as might be seen by specializing the discussion of Sect. 4.1. Alternatively, 
the paths may be viewed as oriented paths on a graph drawn from a square grid in the 
upper-half plane by keeping the horizontal sides and the diagonals of the squares, with 
the two diagonals in each square viewed as passing on top of each other with no vertex at 
their crossing point in the plane. The graph is implicitly endowed with some orientation 
from left to right so that the paths always have strictly increasing abscissas. On such 
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an oriented graph, we may apply the so-caUed LGV (Lindstrom-Gessel-Viennot) lemma 
||2^ , which states that H^ = deto<i,j<n -^i+j enumerates configurations made of n + 1 
paths connecting the set of points (/j)j=:o,-n with coordinates (— j, 0) to the set of points 
(Oj)j=o,---,n with coordinates (j + 1,0) with no intersections at vertices of the graph. In 
addition to the above weights for up-, down- and level-steps, each configuration receives 
a ±1 factor equal to the signature of the permutation a of {0, ■ ■ ■ ,n} characterizing 
its connections (namely Ij is connected to Oa-(j))- Now it is easy to see that, for such 
configurations of non-intersecting paths, the part of the paths outside the central strip 
is entirely fixed to be a set of straight lines going up on the left of the strip and straight 
lines going down on its right, contributing an overall factor /^"^('^+i)/2 (^ggg pjg ]^q^ 
The only freedom comes from the possibility of crossings (in the plane) along the two 
diagonals of a given square in the strip (since these diagonals do not intersect on the 
graph). Moreover, two such crossings cannot take place on two neighboring squares 
in the strip. An acceptable crossing configuration may therefore be viewed as a hard 
dimer configuration on the segment [0, n], as shown in Fig. 10. Each dimer receives a 
weight W = —{—gsR^''^)^ = —g^R^ with a minus sign for each crossing, equal to the 
signature of a transposition component in the permutation a. In conclusion, the Hankel 
determinant Hn is simply given by 

7^(T^ + l) 
Htt^ ^ R 2 Zhard dimora (6.10) 

on [0,n] 

where Zhard dimers Is the generating function for hard dimers on the segment [0, n] with 

on [0,Ti] 

weight W per dimer. Now it is a classical result that 

^ hard dimers ^ — ^ -^j ( D. 11) 

r„T /-li^.ArT-l-l-l^. V / 

on [U.riJ 



1 


1 


— 


yn+2 


(1 + y) 

1 


n+l 


1 


-y 


y + i 


+ 2 ' 





with the parametrization 

W = J . (6.12) 

y+'y + 2 

Since W — —g^ R^, y is precisely the solution of (|6.7|). 

As for the Hankel minor H^, it is easily identified in the hard dimer language as 
(see Fig. 11 for an illustration) 

~ n(n + l) „ ,„ n'^+n + 1 

Hn ^ [n -\- 1)»S -ft 2 Zhard dimora — gsR R ^ Z tiard dimors . (6.13) 

on [0,n] on [0,n, — 1] 



Expressions ( |678| ) follow immediately via ( pr2| ) and (|1.3| ). 

6.2. Quadrangulations 

Let us now come to the case of quadrangulations, characterized by gk = gA^k,4,- We 
have 5" = in this case so that only paths of even length contribute to P or P^ . The 
only non- vanishing Aq coefficients are Aq and A2, from (|1.5|) , given by 



Ao = R-2g4R^ = l + g4R^, A2 = -g4R (6.14) 
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Fig. 11: Correspondence between configurations of non-intersecting paths 
on the graph of Fig. 9 where the last exit point is shifted by one unit to 
the right, and hard dimers on a segment. We distinguish between the cases 
where the uppermost point (1, n -f 1) is attained (bottom) or not (top). 



with R now given, from ( p..6|) , by 



i? = 1 + 3^4 ^ ■ 



(6.15) 



As quadr angulations belong to the class of maps with even face degrees, we use the 
results of Sect. 5. The non- vanishing S^'s are 



So = 1 — (74 -R , Bi = B-i = —g4^ R 
so that the characteristic equation reads 

1' 



:i-g4R')-9^R'[y + 



y 







By ( |5.11| ), the two-point function reads 

(l-y-)(l-y-+3) 



Rr, — R 



(l-|/"+i)(l-y'^+2) 
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(6.16) 



(6.17) 



(6.18) 






\+9,R 




-21 -4-2 2 4 6 



2j+2 



Fig. 12: Interpretation of formula (6.20) for F2i+2j in the case of quad- 
rangulations. The weights per step are different (and allow for a level- step) 
in the central strip as shown. 
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2n+2 



+ 

Si 
<N 



Fig. 13: Correspondence between configurations of non-intersecting paths 
on the graph of Fig. 12 and hard dimers on a segment of length 2?! + 1. 

which matches Eq. (4.10) of [^ (up to the change of notation y -^ x and n — )■ n + 1). 
Again, we may give a combinatorial hard-dimer interpretation to these results as 
follows. The generating function Fn vanishes for n odd and we have 



2n 



Ao P+(2n; R, 0) + A2 P+(2n + 2; R, 0) = Aq P^cat^ + A2 R'^+'catn+i (6.19) 



where cat^ = (^)/(^ + 1) are the Catalan numbers. In particular F2i+2j may be 
interpreted as enumerating paths from, say (— 2z, 0) to (2j -|- 2, 0) on the graph made 
of the restriction in the upper half-plane made of a tilted square grid, completed by 
horizontal segments in the central strip between abscissas and 2 (see Fig. 12). Again 
the graph is implicitly oriented from left to right so that paths have increasing abscissas. 
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Each up- or down-step of the path receives a factor ^/R except those in the central 
strip which receive instead a weight yR\/A2 = iy^-R, while the horizontal paths 

receive a weight Aq = 1 -\- (74 -R^. Using again the LGV lemma p^, the quantity hn = 
deto<ij<n F2i+2j enumerates sets of n-|-l non-intersecting paths from points (/j)j=o,-n 
with coordinates (— 2j, 0) to points (Oj)j=o,-ri with coordinates {2j + 2, 0). Again these 
paths have fixed ascending and descending parts on both sides of the central strip, 
contributing an overall factor R^(^+^) and the only freedom comes from the central strip 
where the j'-th path (numbered to n from bottom to top) connects (0, 2j) to (2, 2j) 
either via a horizontal step with weight 1 -t- g^B?, or by a two-step sequence passing 
either by (1, 2j - 1) (if j > 1) or by (1, 2j - 1), with a total weight {ly/glRY = -g^R^- 
Since the paths are non-intersecting, these (up or down) two-step sequences cannot be 
adjacent in the strip and their vertical positions define a hard dimer configuration in 
[0, 2n + 1] (see Fig. 13). Extracting an overall factor (1 + g^R'^)'^'^^ , each dimer receives 
a weight 



W 



g^R' 



i + g4R^ 



so that we have 



on [0,2n + l] 



(6.20) 



(6.21) 



K 



2i+2j+2 




Fig. 14: Interpretation of formula (6.20) for ^2^-1-2^+2 in the case of quad- 
rangulations. Note that the central strip is shifted by one unit with respect 
to Fig. 12 

As for hn = deto<i,j<n -^2i+2j+2, it enumerates sets of n-|-l non-intersecting paths 
from points (/j)j=o,-n with coordinates (— 2j, 0) to points (Oj)j=o,-n with coordinates 
(2j + 4, 0) with now a strip between abscissas 1 and 3 (see Fig. 14), easily transformed 
into a hard dimer configuration in [0, 2n -|- 2] (see Fig. 15). We deduce from this picture 
the formula 

/i(i) = i?("+i)'(l + £,4i?2)'^+i^,_^ ^;„_^ . (6.22) 

on [0,2n + 2] 
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Fig. 15: Correspondence between configurations of non-intersecting paths 
on the graph of Fig. 14 and hard dimers on a segment of length 2n + 2. 



Using the explicit form (|6.11|) , this leads eventually to 



-K2n+1 — TTrT' '"•- ~ ^ 



,(0) 



"'n-1 "-n-l 



R 



h 



(0) 1,(1) 



h. 



2n 



' ,_(1 



n-1 



R 



'I _ y2n+2\(l _ y2n+3-) 



(6.23) 



h^n\ hl% "(l-y2n+l)(l_^2n+2) 



where y is precisely the solution of ( |6.17|) . The above equations are summarized into 

(iH). 



7. Conclusion and discussion 

To conclude, we discuss further connections of our results raising some open ques- 
tions. Sect. 7.1 is devoted to the manifestations in our context of the intimate relation 
between continued fractions and orthogonal polynomials. Sect. 7.2 makes the connection 
with the matrix integral approach. 

7.1. Convergents and orthogonal polynomials 

There is a deep connection between continued fractions and orthogonal polynomials. 
Let us indeed consider the family of polynomials (Qn(^))n>o defined by 



Qn{z) 



1 



det 



Hn—l 0<i,j<n 



F^+j]z^ 



(7.1) 



where {oij ] Oj) denotes the matrix with elements Oj in the last row and elements Oij 
in all previous rows. It is a classical exercise to check that qn{z) is a monic polynomial 



36 



of degree n, that the family is orthogonal with respect to the scalar product defined by 
{z'^jZ'^) = Fn+rm with (q'n(^) 7 Q'm(^)) = Sn,mHn/ H^-i^ and that it satisfies Favard's 
three-term recurrence 

Zqniz) = Qn+liz) + Snqniz) + RnQn-liz) (7.2) 

with the initial data qo{z) = 1, q-i{z) = 0. The connection with continued fractions 
comes from the fact that the reciprocal of q-niz) appears as the denominator of the nth 
convergent of the J- fraction ( |1 . 1| ) , namely 



Pn{z) _ 1 

qn{z) RlZ^' 
i — ^0^ 



(7.3) 



B~z^ 

1 - bn-lZ- 



1 - Sr,.Z 



where qn{z) = z^qnil/z) and where Pn{z) are the so-called numerator polynomials. 
•Pniz) and qn{z) also appear within an expression for the nth truncation (|2.2|), see for 
instance [|T|,|^. 

In the context of maps, the dth convergent of (|1 . 1|) may be interpreted as the 
generating functions for rooted maps where every vertex incident to the root face is 
at a distance lesser than or equal to d from the origin of the root edge. As for the 
dth truncation (p.2| ), its interpretation is given in Sect. 3.3 via (|3.13|) as a generating 



function for pointed rooted maps with a control both on the root degree and on the 
distance from the origin to the root face. Truncations are known to admit a simple 
expression in the case of quadrangulations |2^] , namely 



1 l-{W-l)y\ 



.yd+l 



_W ^ ^^-^ " (7 4) 

-. _ Rd+2Z 
~ 1 

with K and y as in Sect. 6.2 and W = {\ — Vl — ^Kz^^ j {2Rz^\ We may wonder if this 
expression generalizes to arbitrary (/^'s. We provide in Appendix C a general expression 
for the orthogonal polynomials qn{z\ Unfortunately we lack an expression for the 'Pn{z\ 
as they do not admit a formula similar to ( |7.1| ) . 

1.2. Connection with matrix integrals 

Another classical and fruitful approach to map enumeration problems is via matrix 
integrals. Let us now comment informally on its connection with our present results. 
Matrix integrals give a simple expression for the all-genus generating function 

F4N] = J2n~^^FJ^''^ (7.5) 

h>0 
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where Fn denotes the generating function for rooted maps of genus h with a root face 
of degree n, with face weights {gk)k>i as before. In particular, we have Fn = Fn = 
hniAr^oo Fn [N] . We have the matrix integral representation 

^,,,,_ J dMTr{M^)eM-NTrV{M)) 
"^^ N J dM expi-NTtViM)) ^ ' 

where dM denotes the Lebesgue (translation-invariant) measure over the space oi NxN 
hermitian matrices and 

9 h 

nx)^y-E^^T- (7-7) 

k>l 

The terms gux^ act as a "perturbation" of the quadratic potential x^/2 corresponding 
to the well-known Gaussian Unitary Ensemble. In all rigor, expression (|7.6|) must be 
understood as a power series in the g^s and A^, whose precise definition is beyond the 
scope of this section. 

We now briefly discuss the usual approaches for studying (|7.6| ), beside the loop 
equations already mentioned in Sect. 3.2. The original "physical" approach is the so- 
called saddle-point or steepest descent method 0. It consists in remarking that the 
integrands in ( [7.6|) only depend on the eigenvalues of M, and observing that for large 
A^ the dominant contribution comes the "equilibrium" continuous distribution of eigen- 
values. It yields 

Fn= j d\\^p{X) (7.8) 



where p(A) denotes the density of eigenvalues. Therefore, in the saddle-point picture 
our moments Fn are precisely those associated with the spectral measure. 

This has to be contrasted with the usual method of orthogonal polynomials in 
random matrix theory Pj. There, we consider the family of polynomials {ql (A))i>o, 

with ql (A) monic of degree z, orthogonal with respect to the scalar product deflned by 
(A"^, A"^) — J dX A'^+"^ exp(— A^ V{X)). In other words, these orthogonal polynomials are 
deflned with respect to the "A^ = 1" eigenvalue density while those of ( [7.1| ) are deflned 
with respect to the "A^ — )■ oo" eigenvalue density. Favard's theorem states that we still 
have a three-term recurrence 

A<,r'(A) = 4:1 W + Sr> ,f '(A) +iii~' <,._,(A) (7.9) 

where {S^ )i>o and {Ri )i>i (with Rq =0) are associated with the scalar product 
at hand, for instance they may be expressed via Hankel determinants of the moments of 
the measure dX exp{—NV{X)). Now, because of the speciflc form of the scalar product 
(■,■), we also have the relations f^ 



«. ' = 1^+E9'C^(*-1) 



A^ 

'^-^ (7.10) 

k>i 
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Here Z^-\k) denotes the generating function of three-step paths from (0,i) to (/c,j) 

with a weight Rm per down-step (t, m) — )■ (t -|- 1, m — 1) and a weight Sm per level- 
step (t, m) — )■ {t + l,m). Remarkably enough, these equations look very similar to (|2.4|) , 

and allow to identify R^ and S^ with generating functions for mobiles (respectively 
half-mobiles) as defined in this paper with, however, a somewhat mysterious and non- 
conventional weight [m/N) for labeled vertices with label m. In terms of orthogonal 
polynomials, Fn[N] is given by 



N^ ^fHA),,f)(A)^ ~N 






i=0 (/i lAj,y. yA) i=o 



We may further connect the relations ( 7.10 ) to known mobile generating functions 



by considering the limit N ^ oo. In this limit, the quantity u = i/N is treated as 
a continuous variable, and R\ and Sn become at leading order in 1/A^ smooth 
functions R{u) and S{u). Then, in this limit ( |7.1C1| ) yields 

R{u) = u + Y.9kVW)P-i{k - 1; R{u), S{u)) 

"-" (7.12) 

S{u) = Y,9kP{k~l-R{u),S{u)) 

fc>i 

involving paths which receive homogeneous weights R{u) and S{u), irrespectively of the 
ordinates. Setting w = 1, we recover precisely (|2.5|). For general w, R{u) and S{u) are 
generating functions for mobiles and half-mobiles with an extra weight u per labeled 
vertex, as encountered in Sect. 3.1 and Appendix B. In the limit A^ — > oo, the sum in 
( [7.11| ) becomes an integral, leading to the expression 

1 

duP{n;R{u),S{u)) (7.13) 



which is consistent with the discussion of Sect. 3.1. 

For finite A^, a natural question is whether the unconventional weight of labeled 
vertices, as well as the form (|7.11|) of the all-genus generating function may be given a 
direct combinatorial interpretation. This could open the way to get explicit formulas 
for discrete distance-dependent two-point functions in maps of higher genus. 

Acknowledgments: We thank C. Krattenthaler for pointing out the connection with 
Schur functions. Part of this work was completed at the Centre Emile Borel of Institut 
Henri Poincare. 
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m ■^^ I 
1=1 



Fig. 16: The decomposition (a) of a pointed rooted map into slices by 
cutting it along all leftmost geodesic paths emerging from vertices incident 
to the root face. Slices are of two types: those (b) with a basal edge of type 
{m,m), counted by 5"^, and those (c) with a basal edge of type (m,m — 1), 
counted by Rm- Note that edges of type (m — 1, m) delimit empty slices. 
The actual depth of a slice may be less than m since its two boundaries 
may merge before the origin {si and r^ counting the slices with depth i). 



Appendix A. Decomposition of maps into slices 

So far we used mobiles as the natural framework to introduce continued fractions 
in map enumeration problems. Still, the explicit recourse to mobiles is not stricto sensu 
necessary to prove the basic continued fraction expansion (|1 . 1| ) and more generally the 
interpretation ( ^.131 ) of truncations ( |2.2D , when the weights Rm and Sm have their 
original interpretation of Sect. 1 as map generating functions. This relation can be 
understood alternatively as the result of some appropriate decomposition of maps which 
can be defined without reference to the underlying mobile structure even though, in 
practice, the two are intimately related. 

We start again with pointed rooted maps enumerated by fn;d, i-e. maps with a root 
face of length n, with the origin of the root edge at distance d from the origin of the 
map and with all the other vertices incident to the root face at a distance larger than 
or equal to d. As before, we choose the root face for external face and the clockwise 
sequence of distances of its incident vertices defines a three-step path of length n with 
endpoints at ordinate d, and which stays above d. For each vertex incident to the root 
face, we draw the leftmost geodesic path from this vertex to the origin of the map. 
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Cutting along these paths decomposes the map into a number of connected domains 
that we caU slices (see Fig. 16). These shces are maps whose boundary is made of 
one basal edge originaUy incident to the root face and of two leftmost geodesic paths 
merging at some apex (possibly different from the former origin of the map). Note 
that the leftmost geodesies may themselves contain edges originally incident to the root 
face. The basal edge of a non-empty slice is an edge originally of type (to, m) or of type 
(to, to — 1) (counterclockwise). Indeed, the edges of type (to — 1, to) give rive to empty 
slices as they lie on the leftmost geodesic path emerging from their endpoint at distance 
m. Non-empty slices are therefore associated only with the down- and level-steps of 
the three-step path above. Since two consecutive leftmost geodesies may merge before 
reaching the origin of the map, the actual length of the leftmost geodesies is reduced 
in the slice by some value m — z corresponding to the length of their common part. 
This leads us to associate a weight Sm = Xl^Io ^^ (respectively Rm = YllLi ^j) ^o each 
level- (respectively down-) step starting at ordinate to in the three-step path, where Si 
(respectively r^) are the generating functions for slices of depth i, i.e. with a boundary 
made of a basal edge and two leftmost geodesies of the same length i (respectively of 
length i — 1 and i) (see Fig. 16). Here again, we have to add to ri, hence to Rm, a 
conventional factor 1 which accounts for the case where an edge of type {m,m — 1) 
incident to the root face would be the boundary of an empty slice, which happens when 
its endpoint at distance to — 1 is the only vertex at distance to — 1 from the origin 
adjacent to its endpoint at distance to. On the contrary, each edge of type (to, to) 
incident to the root face is the basal edge of a non-empty slice. We therefore recover 
precisely the generating function ZJ~^(n) defined in Sect. 2.1, which, from the above 
decomposition procedure, enumerates all maps in /^.(^, as well as all maps in /„.j for 
i < d, since the concatenation of slices produces maps whose origin may remain at a 
distance less than d. We end up with the desired relation Fn;d = ^ddi''^)^ ^^^ ^^^ with 
a different interpretation for the weights 5"^ and Rm as slice generating functions. It 
remains to show that these new definitions match their former definition of Sect. 1.2 in 
terms of map generating functions. 

It is easily seen (Fig. 17) that, starting from a slice counted by r^, completing the 
leftmost geodesic of length i — 1 by the basal edge (of type (i — 1,«)) and gluing it 
with the other leftmost geodesic of length i produces bijectively a pointed map with 
a marked edge of type {i,i — 1) with respect to the origin of the map. This marked 
edge is promoted as root edge. Note that, after gluing, we are apparently left with 
a marked geodesic but it is by construction the leftmost geodesic starting with the 
root edge (oriented from i to i — 1) and ending at the origin. It is therefore uniquely 
determined and may be erased without loss of information. This bijection guarantees 
that our new definition of r^ (and Rm) matches the definition of Sect. 1. As for Sn, 
we see (Fig. 17) that, upon gluing two slices in Sn (respectively in Si and Sj for some 
i < n and j < n), we get a pointed map with a marked edge of type (to, to) for some 
m < n (with to = max(i, j)). This marked edge is oriented with, say, the first slice on 
its right, and promoted as root edge. Again, we have after gluing two marked leftmost 
geodesies (starting respectively with the root edge and the reversed root edge) which 
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Fig. 17: Gluing two slices in Si and Sj and identifying pairwise boundary 
edges creates a pointed rooted map of type (m — )■ m) with m = max(z, j), 
implying S"^ = T„ with T„ defined as in Sect. 1. Identifying pairwise 
boundary edges of a slice in r^ creates a pointed rooted map of type {i — )■ 

may be erased without loss of information. Hence 5"^ = Tn with Tn defined as in Sect. 1. 



Appendix B. Proof of ( pTSl) 

Let us start by showing the following identity, valid for any n > and k >2 

g=0 



^(RY,P^in + q;R,S)P{k-2-q;R,S)]=Pin;R,S)^[VRP-i{k-l-R,S)) 



+ VRP-i{n-R,S)^P{k-l-R,S) 

(B.l) 
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At this stage R and S are arbitrary coefficients. Attaching a weight x'^y'^ ^ and summing 
over n > and k > 2, the r.h.s produces the combination 

d f 7V-i{y) \ d f TT{y)-l \ 

"^") dR [-^) + "-^^"^ dR [-^) ' ^^-'^ 

where we introduced the generating functions 

1 



n{z) = ^P(?i;i?,5); 



n>0 \/^ 

1 — Sz — ^/k{z) 

V JTL JT — iyn, JTL, o )^ ^Z 

n>0 



7r_i(z) = J]/RP_i(n;i2,5) 



(B.3) 



2 = ;2- 

2Rz'^^l^) 



with ^(2) as in ( |3.8| ). The last formula comes from the identification 7r_i(2:) 
27r"'"(2)7r(2), where 



Y.P^in;R,Sy,^ '-%-f^^ 



n>0 



As for the l.h.s in (|B.1|) , setting s = n + q and t = k — 2 — q, we have to compute 

fc-2 

E E E ^^(^ + ^; ^' ^)^(^ - 2 - ?; i?, ^)a;-y'=-2 

n>0 fc>2 (j=0 

s 

s>0t>0g=0 (B.5) 

= E E ^ 2/*i^+(s; i?, S)Pit; R, S) 

s>0 t>0 ^ 

xn+ix) -yn+{y) 

= 7t{y) . 

x-y 

Proving (|B . 1|) therefore reduces to checking the relation 

d f^ x7r+{x)-y7r+{y) \ d f n-i{y) \ d fr:{y)-l 



dR \ X — y J dR \ y J OR \ y 

(B.6) 
which, knowing the explicit forms of all the involved generating functions, is a straight- 
forward task. 

Multiplying ( [B.l| ) by 5k,2 — 9k, summing over k >2 and exchanging the sums over 
k and q leads to 

— J2A,P^in + q;R:S) = Pin;R,S)—u + VRP-iin;R,S)—v (B.7) 

q=0 
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where we introduced 

u = R-^gkVRP-i{k-l;R, S) 

"-" (B.8) 

v = S-Y,9kP{.k-l-R,S) 

k>i 



and where Aq is defined as in ( p. . 5| ) . Note that we added for convenience a trivial 
constant term gi in the definition of v, which disappears after derivation with respect 
to R. 

Similarly, one can easily prove 

— iRj2p+{n + q;R,S)Pik-2-q;R,S)] = P{n;R, S)— (VRP-iik - 1; R, S)^ 

+ VRP-i{n;R,S)^P{k~l;R,S) 

(B.9) 
by checking the identity 

d f X7r+{x)-yn+{y) \ d fTV-i{y)\ d f 7r{y) - 1 



dS \ X — y J dS \ y J dS \ y 

(B.IO) 
This implies the relation 

— Y, AqP+{n + q- R, S) = P{n- R, S) —u + ^^P-i{n- R, S) —v (B.ll) 

q=0 

with u and v as above. Equations ( [B.?] ) and ( [B.ll| ) may be rewritten in terms of 
differential forms 

dl^AqP+{n + q]R,S)] =P{n;R,S)du + VRP-i{n;R,S)dv . (B.12) 

When V = 0, Eqs.( p^) match precisely the equations satisfied by the mobile generating 
functions R{u) and S{u) with a weight u per labeled vertex. Therefore, by restricting 
( [B.12| ) to the line f = 0, we deduce 

J oo 

— J2 M^)P^i^ + 9; R{u),S{u)) = P{n- R{u),S{u)) (B.13) 

'" q=0 



with Aq{u) defined as in ( p..5| ) with R, S replaced by R{u),S{u). This establishes 
Eq. (|3.3|) upon integrating with respect to u and noting that both sides of ( p.3|) vanish 
at u = {Aq{0) = since -R(O) = and Fn{0) = from its map definition). 
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Appendix C. An expression for the orthogonal polynomials 

Using the notations of Sect. 4, we have qn{z) = — det I ^ ) where z^ is the 

Hn-l \ Zn 

row vector (^"')q< < • Inspired by ( [4.2| ), we may write 



i; Un-T- 



(C.l) 



By Proposition 11 of fl^, the inverse matrix of T = {P^{j; R, S))ij>o is given 
by the matrix of coefficients of orthogonal polynomials associated with the contin- 
ued fraction (|2.3|). The jth such polynomial is Uj{{z — S)R~^^'^), where Uj{z) de- 
notes the jth Chebyshev polynomial of the second kind (defined for instance via 
Uj{2cose) sin^ = sin(i + 1)^). Hence we have z„ ■ T"! = {Uj{{z - S)R-^/^))^^^^^. 
Passing to determinants in (|C.1|) , expanding the middle r.h.s determinant over the last 
row and using the known expressions for det T^ and Hn , we obtain 

(C.2) 



m=0 



m=0 ^ 



"J \ ^P2p\'^P,'n 5 ^ 



VR J SP2p(Ap,n,x) 

where Ap'^ = 'nP~^m denotes the partition with p — 1 parts of size n and one part of 
size m. The h- or Weyl formulas allow to further reexpress qn{z) in terms of p x p 
determinants, where only one row depends on z. 
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